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Abstract 

A transformation is derived which takes the Lorenz integrable system into the well- 
known Euler equations of a torque-free rigid body about a fixed point, i.e. the famous 
motion a la Poinsot. The proof is based on Lie group analysis applied to two third- 
order ordinary differential equations admitting the same two-dimensional Lie symmetry 
algebra. Lie's classification of two-dimensional symmetry algebras in the plane is used. 
If the same transformation is applied to the Lorenz system with any values of the 
parameters, then one obtains Euler equations of a rigid body about a fixed point 
subjected to a torsion depending on time and angular velocity. The numerical solution 
of this system yields a three-dimensional picture which looks like a "tornado" the cross- 
section of which has a butterfly-shape. Thus Lorenz's butterfly has been transformed 
into a tornado. 
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1 Introduction 



The motion of a heavy rigid body about a fixed point is one of the most famous problems of 
classical mechanics P]. In 1750 Euler derived the equations of motion, which now bear 
his name, and described what is nowadays known as the Euler-Poinsot case because of the 
geometrical description given by Poinsot about hundred years later It was Jacobi jl] 
who integrated this case by using the elliptic functions which he had developed (along with 
Legendre, Abel and Gauss and mastered [H] - we have translated this fundamental text 
into Italian and commented extensively '7]. 

More than 200 years later, in 1963, a paper was published [Hj in which was presented a 
system of three ordinary differential equations. The author considered a hydrodynamical 
system developed by Rayleigh and reduced it by applying a double Fourier series as in 
[To] . Thus he obtained what nowadays is the famous Lorenz system pi]. Three parameters 
are part of the Lorenz system. For particular values of those parameters the Lorenz system 
can be integrated in closed form by means of Jacobi elliptic functions jT2|. We call this 
system the Lorenz integrable system. 

In January 2001 the first Whiteman prize for notable exposition on the history of mathemat- 
ics was awarded to Thomas Hawkins by the American Mathematical Society. In the citation, 
published in the Notices of AMS 48 416-417 (2001), one reads that Thomas Hawkins ". . . 
has written extensively on the history of Lie groups. In particular he has traced their origins 
to [Lie's] work in the 1870s on differential equations . . . the idee fixe guiding Lie's work was 
the development of a Galois theory of differential equations . . . [Hawkins's book ^Hj] high- 
lights the fascinating interaction of geometry, analysis, mathematical physics, algebra and 
topology . . ." . Also Hawkins had established "the nature and extent of Jacobi's influence 
upon Lie" [H]. 

In the Introduction of his book ^2] Olver wrote that "it is impossible to overestimate the 
importance of Lie's contribution to modern science and mathematics. Nevertheless anyone 
who is already familiar with [it] ... is perhaps surprised to know that its original inspirational 

2 



source was the field of differential equations" . 

Lie's monumental work on transformation groups, jTHj, |IZI and ^Hl, and in particular con- 
tact transformations ^Hl, led him to achieve his goal |2()j . 

Many books have been dedicated to this subject and its generalizations [21]; jUj; 

EH, Ezi, m and m)- 

Lie group analysis is indeed the most powerful tool to find the general solution of ordinary 
differential equations. Any known integration technique can be shown to be a particular 
case of a general integration method based on the derivation of the continuous group of sym- 
metries admitted by the differential equation, i.e. the Lie symmetry algebra. In particular 
Bianchi's theorem ([10], J3]) states that, if an admitted n-dimensional solvable Lie sym- 
metry algebra is found, then the general solution of the corresponding n*''-order system of 
ordinary differential equations can be obtained by quadratures. The admitted Lie symmetry 
algebra can be easily derived by a straightforward although lengthy procedure. As computer 
algebra software becomes widely used, the integration of systems of ordinary differential 
equations by means of Lie group analysis is becoming easier to perform. A major drawback 
of Lie's method is that it is useless when applied to systems of n first-order equations, be- 
cause they admit an infinite number of symmetries, and there is no systematic way to find 
even an one-dimensional Lie symmetry algebra, apart from trivial groups like translations in 
time admitted by autonomous systems. One may try to derive an admitted ra-dimensional 
solvable Lie symmetry algebra by making an ansatz on the form of its generators. 
However, in |^ we have remarked that any system of n first-order equations could be trans- 
formed into an equivalent system where at least one of the equations is of second-order. 
Then the admitted Lie symmetry algebra is no longer infinite-dimensional, and nontrivial 
symmetries of the original system could be retrieved [SI]. This idea has been successfully 
applied in several instances (jSIj, [22]; [22] and [21] )• Also in |23] we have have shown that 
first integrals can be obtained by Lie group analysis even if the system under study does 
not come from a variational problem, i.e., without making use of Noether's theorem jSHj- If 



3 



we consider a system of first-order equations and, by eliminating one of the dependent vari- 
ables, derive an equivalent system which has one equation of second-order, then Lie group 
analysis applied to that equivalent system yields the first integral(s) of the original system 
which do(es) not contain the eliminated dependent variable. Of course this requires that 
such first integrals exist. The procedure should be repeated as many times as there are 
dependent variables in order to find all such first integrals. The first integrals correspond to 
the characteristic curves of determining equations of parabolic type which are constructed 
by the method of Lie group analysis. We remark that interactive (not automatic) programs 
for calculating Lie symmetries such as |^ and |SHI are more appropriate for performing this 
task. 

We have briefiy sketched three apparently unrelated subjects. In this paper we show that the 
Lorenz system and the Euler equations are actually related by means of Lie group analysis. 
In inn] we applied Lie group analysis to a third-order differential equation, which is equiva- 
lent to the Lorenz integrable system, and obtained a two-dimensional Lie symmetry algebra, 
which we then used to integrate the Lorenz integrable system in terms of Jacobi elliptic 
functions. Here we show that the same Lie symmetry algebra is admitted by a third-order 
differential equation which is equivalent to the Euler equations of a torque-free rigid body 
moving about a fixed point. Then a transformation is easily derived by which the Lorenz 
integrable system becomes the Euler equations of a torque-free rigid body moving about a 
fixed point. Thus, it can be stated that "the Lorenz integrable system moves a la Poinsot". 
If the same transformation is applied to the Lorenz system with any value of parameters, 
then one obtains the Euler equations of a rigid body moving about a fixed point and sub- 
jected to a torsion depending on time and angular velocity. The numerical solution of this 
system yields a three-dimensional picture which resembles a "tornado" the cross-section of 
which has a butterfiy-shape. By means of our transformation Lorenz's butterfly becomes a 
tornado. 

In the last section the relationship between Lie group analysis and first integrals [HSl is ex- 
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emplified by considering the Euler equations of a torque-free rigid body about a fixed point 
and the Lorenz integrable system. 

2 Butterflies and tornadoes 

Consider the Lorenz system |H]: 

x' = a{y-x), (1) 

y' = —xz + fx — y, (2) 

z' = xy — bz, (3) 

where a, b and f are parameters (a prime denotes differentiation with respect to r). This 
system can be reduced to a single third-order ordinary differential equation for x jlT] , which 
admits a two-dimensional Lie symmetry algebra if (7 = 1/2, 6=1 and f = 0. System (PO)) 
becomes 

^' = (4) 

y' = -xz-y, (5) 

z' = xy — z. (6) 
The corresponding third-order equation is: 

2xx"' — 2x'x" + 5xx" — 3x'^ + 2x^x' + 3xx' + + a;^ = 0, (7) 
and admits a two-dimensional Lie symmetry algebra L2 with basis: 

Xi = dr, X2 = (^Or - ^Xd,^ . (8) 

A basis of its differential invariants of-order < 2 is given by: 

^=(^x' + I) x-\ ^ = [x" + ^x' + I) x~\ (9) 
Equation ((Zj) is reduced to the following first-order equation: 

(V^-202)^ = -2V^0-0, (10) 



which can be easily integrated |40j to give: 



Cl, 



(11) 



(1 + 2^)2 

where ci is an arbitrary constant. Substitution of x and its derivatives into (jllll yields a 
second-order ordinary differential equation 

1 + 4 [x" + fx' + x'^ -A{x' + \x^ ^ x~'^ 



(x^ + 2x" + 3x' + xf x-^ 



Cl, 



(12) 



which admits the Lie symmetry algebra L2. Lie's classification of two-dimensional algebras 
into four canonical types [201 allows us to integrate (fT^ by quadrature if we introduce the 
canonical variables: 



-r/2 



U 



X 



(13) 



which transform equation ()12|) into 



4m- 



d^u 
df2 



2n— -4 
df 2 



d-u 
dt; 



- 1 



Cl, 



(14) 



and operators (jHI) into 



Xi = (9^, X2 = vdy + udu . 
Then the general solution of (fT^ can be easily derived [20] to be: 



-Cl =F 202^^ — C2'U'' 



-1/2 



du 



+ C3, 



(15) 



(16) 



with C2 and C3 arbitrary constants. This solution which involves an elliptic integral has 
already been obtained by Sen and Tabor |3T] by means of a lengthier analysis. 
The Euler equations describing the motion of a heavy rigid body about a fixed point with 
no torsion are 

{B-C) 



V 

<i 



A 

jC-A) 
B 

C 



qr, 
pr, 
pq, 



(17) 
(18) 
(19) 



with A, B and C being the principal moments of inertia, and q{t) and r{t) the compo- 
nents of the angular velocity (a dot denotes differentiation with respect to t). This system 
can be reduced to a single third-order ordinary differential equation for, say, p, viz 

d?p _ dpd?p _ A{C-A){A-B) 3dp _ 
^dt3 dtdt^ EC ^dt"^' ^^"^ 

which admits a two-dimensional Lie symmetry algebra C2 with basis: 

Ti = dt, T2 = tdt - pdp . (21) 

The two Lie symmetry algebras L2 and C2 that we have found are actually the same, i.e. 
Type IV in Lie's classification [20] ■ Therefore they are linked by a transformation which 
takes (t, x) into {t,p). Prolongation to the second-order of the two equivalent Lie symmetry 
algebras yields a transformation which takes the system (fT7j) - ()19j) into the system (^-(jB)) as 



log 

pt 



2 



C-B 2 
-qrt , 



2A 
C-B 



B C ^ 



2A 

with the following condition on the momenta of inertia 

{A-B){A-C) 1 



(22) 



(23) 



BC 4 

A slightly more general condition could have been considered if one replaces 1/4 with k/4 {k 
an arbitrary parameter). If one derives B from (I^Hj) by assuming A — C>0 and AA — SC > 0, 
i.e. 

_ 4A{A - C) 

^ - AA-3C ' ^^^^ 
then the transformation (j22j) turns into the following 
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pt 

^ {2A-C)(2A-^C) 
^ 2A(4A-3C) ^ ' 

(2A - C){2A - 3C) [4A2g2 _ (4^ _ ac)^^] ^ 
^ ~ 8^2 (4A - 3C)2 ^ ' 

and the system (fT7j) - (fT^ assumes the form 



(25) 



(2A-C)(2A-3C) , , 

3C - 

9 = — — Vr. (27) 

(28) 



4^-3^ 



We also derive the inverse transformation, i.e. 



t = 2e-^/2^ 
p = xe"^/^, 

2Ji2A - C)i2A - 3C)(v/^^T^ +z)' 



Ae^/^J^y' + + 



' {2A - C){2A ~ ?>C) 

which takes the system (jl))-® into the system (fT7j) - (fT9|) after substituting B as in plj) . 
If one apphes the transformation ()29|) to the general Lorenz system (^Q)-©, then the following 
equations are obtained 

2{2A-C){2A-'iC)a , ,P 

3C - 4A ,~ ^ , 4A2g2 _ UA - 3C)2r2 q 
Q = nr + ib-l) 

, . 2(4A - SCy^A pr- 

{2A - C) {2A - 3C) [AA^q^ + {AA - 3Cyr^] ' ^ ' 

. _ A ~ 4AV - (4A - 3C) V r 

~ 4A - 3C^^ ~ ^ ~ MA2g2 + (4A _ 3C)2r2 t 

, . 8(4A - 3C)A^ p_q 

{2A - C) i2A - 3C) [AA^q^ + {4A - 3Cyr^] ' ^ ^ 
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They can be interpreted as the Euler equations of a rigid body moving about a fixed point 
and subjected to a torsion which depends on time t and angular velocity [p, q, r) in the 
body-frame reference. Also the momenta of inertia are linked by relation (jSl}. 
To our knowledge such a system has never been described. 

If we use Maple V in order to draw a three-dimensional plot of system (j3(J|) - (j3!2|) . then a 
"tornado", the cross-section of which resembles a butterfiy, is obtained (FIG. 1). The usual 
values for the Lorenz parameters, a = 10, b = 8/3, f = 28, are imposed. Also we assume 
A = 2 and C = 1. We consider t as it varies in the interval [2, 0.015], which corresponds to 
r G [0, 9.8] approximately. The step size used is 0.00005. 

The butterfly-shape curve is better seen in FIG. 2, which shows the two-dimensional plot of 
p versus q. 

A clearer view of the "tornado" is given in FIG. 3, which shows the two-dimensional plot of 
r versus q. Another view can be found in FIG. 4 which shows the two-dimensional plot of r 
versus p. 

In FIG. 5 and 6 the plot of p versus t is given for two different but close initial values (1 and 
1.01). For relatively large values of t, say t G [2,0.5], the solutions are the same (FIG. 5). 
For small t, say t G [0.05,0.015], a dramatic difference appears (FIG. 6). 
In 1121 it found that one could explicitly derive a first integral of the Lorenz system (jT))- 
(jni in six different instances which correspond to particular values of the Lorenz parameters. 
It is a trivial task to apply those findings to system ()30 p - (jn^ by using the transformation 
fl22|) . The following is the list of the six cases with a relative first integral F for the system 

(1) b = 2a, r arbitrary 



F 




(2)6 



0, a = |, r arbitrary 



F 
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{AA\' - (4A - 3C)V2)(2A - C){2A - 3C) ^ i2A - C)\2A ~ 3Cf 



8A2(4A-3C)2 ^ 3A2(4A-3C) 

(3) 6 = 1, f = 0, arbitrary 

{AA^q^ + (4A - 3Cyry{2A - Cf{2A - SC^ 
~ 4A^{4A - 3C)4 ■ 

(4) 6 = 4, a = 1, f arbitrary 

^ . 4r^ (i6.V- . 32. ^^^^^^ - -l';^;'!^^ '^^'^ - - 

^ (4A2g2 - (4A - 3C)V^)(2A - C)(2A - 3C) , , ^ J2A-C){2A-3C) 
^ A\4A-3Cy ^ * A(4A - 3C) ^^'^^ 

,^ J2A-C)^(2A-3C)^ 2 ,, (4A2g2 _ (4A - 3Cyr')i2A - C)i2A - 3C) 

^ A\AA-3CY ^ A\AA-3CY 

(5) 6 = 1, 0" = 1, f arbitrary 

4 2^ (4A2g2 + (4^_3C)V)2(2A-C)2(2A-3C)2 

= T) r H . 

t^^ AA\AA - 3CY 

(6) 6 = 6ct - 2, f = 2a - 1 



p ^ ^4.^2(1-4.) / P't' - (4A - 3Cfr^){2A - C){2A - 3C) 



64a 32yl2(4A-3C)2 
{2a - 1)2^2 + {2A-C){2A-3C){2~a-l) ^^^^ ^ {2A - C)\2A - 3Cf~a 



4a 2A{AA-3C) 4A2(4A - 3^) 

3 Lie group analysis and first integrals 



In jSo] we showed the apphcation of Lie group analysis in order to obtain first integrals with 
at least one missing variable. Consider Euler equations (fT7j) - (fT!H) . There exist two well- 
known first integrals, i.e. the conservation of kinetic energy and the conservation of angular 
momentum: 

V + Bc^ + Cr'^ = h, (33) 
AV + 5V + CV = J2. (34) 
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We derive p from (fT^ (the method works independently of the chosen equation) as: 

_ Cf 
^~ {A-B)q 

in order to obtain the following two differential equations in q and r, one of first-order and 
one second-order, respectively: 

C{CA)rr 



B{A-B)q 

C{C-A) .2 , ^B-C){A-B) 2 



(35) 



B{A-B)q^''^ AC ^^'^ 

When Lie group analysis of the system (jH^jl - ljHBj) is performed, a linear partial differential 
equation of parabolic structure is generated. Its characteristic curve is given by 

B{A- B)q^ + CiA-Cy 

which is a combination of the two first integrals (j^Hj) and (jMj) . Consequently we introduce 
the new dependent variable s such that 



Q 

in order to obtain the following system 



C{C - A)r^ + s 
\ B{A - B) 



(37) 



s = 0, (38) 

.. ^ C{C~A)tP {C-A){B-Cy {B-C)sr 

^ C{C-Ay + s AB ABC ' ^ ^ 

Equation ()39|) admits an eight-dimensional Lie symmetry algebra (i.e., it is linearizable) if 
either A = C or B = C, i.e. the case of the torque- free Lagrange top (uniform precession). 
If either g or r, one at a time, is eliminated from system (fT7j) - (fT^ . then a similar result is 
obtained, i.e. the other two combinations of the first integrals (jH^ and (jMj) . Indeed the 
elimination of q yields the first integral 

A{A - By + C{B - C)r\ 
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and the elimination of r yields the first integral 

B{B -C)q^ + A{A-C)p^. 

If the same method is applied to the Lorenz integrable system (0))-©, then the following 
two first integrals are obtained 

{y' + z') (40) 

{z-x'^)e\ (41) 

Those first integrals were found in [T? by using the Painleve analysis. 
From (jH)) we have 

z' + z 
X = 

y 

in order to obtain the following two differential equations in y and z, one of first-order and 
one second-order, respectively: 

= -yl±J!l±il, (42) 



y 

- 3y^z - 5y^z' - 2z^ - Az^z' - 2zz'^ 



z" = . (43) 



When the Lie group analysis of the system ()42|) - (j4H|) is performed, a linear partial differential 
equation of parabolic structure is generated. Its characteristic curve is given by 

y +z 

Consequently we introduce the new dependent variable Y such that 



y = VY -z^ (44) 

in order to obtain the following system 

Y' = -2Y, (45) 
Y^ -2Yz^ -3Yz-5Yz' + z^ + z^ + z^z' -2zz'^ , , 

^ = W^) ■ 

The first equation can be easily integrated to yield the first integral (jlO)). 

If we eliminate y from system (jH)-®, a similar procedure provides the first integral PT| . 
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Plot: figrigid01.eps 



Figure 1: The 3-dim plot of the "tornado". 



Plot: figrigid05.eps 



Figure 2: The 2-dim plot of p vs. q. 



Plot: figrigid06.eps 



Figure 3: The 2-dim plot of r vs. q. 



Plot: rigids04.eps 



Figure 4: The 2-dim plot of r vs. p. 



Plot: rigids03.eps 



Figure 5: Varying the initial condition of p in i = 2 by 0.01. The two plots ofp{t),t e [2, 0.5] 
are indistinguishable. 



Plot: fig6.eps 



Figure 6: Varying the initial condition of p in t = 2 by 0.01. There are two different plots of 
p{t),t e [0.05, 0.015]: the darker hne represents the plot of p{t) with initial condition p —1 
in i = 2, and the lighter line the plot of pit) with initial condition p — 1.01 int — 2. 
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